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Abstract In this paper, we propose a variational multiphase image segmen¬ 
tation model based on fuzzy membership functions and L 1-norm fidelity. Then 
we apply the alternating direction method of multipliers to solve an equivalent 
problem. All the subproblems can be solved efficiently. Specifically, we propose 
a fast method to calculate the fuzzy median. Experimental results and com¬ 
parisons show that the L 1-norm based method is more robust to outliers such 
as impulse noise and keeps better contrast than its L2-norm counterpart. The¬ 
oretically, we prove the existence of the minimizer and analyze the convergence 
of the algorithm. 

Keywords Image segmentation; fuzzy membership function; L 1-norm; 
ADMM; segmentation accuracy. 


1 Introduction 

As a fundamental step in image processing, image segmentation partitions an 
image into several disjoint regions such that pixels in the same region share 
some uniform characteristics such as intensity, color, and texture. During the 
last two decades, image segmentation methods using variational methods and 
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partial differential equations are very popular due to their flexibility in problem 
modeling and algorithm design. There are two key ingredients of variational 
segmentation models. One is how to describe the regions or boundaries between 
these regions, and the other is how to model the noise and describe the uniform 
characteristics of each region. 

The Mumford-Shah model a well-known variational segmentation 

model, penalizes the combination of the total length of the segmentation 
boundaries and the L2-norm error of approximating the observed image with 
an unknown piecewise smooth approximation. In other words, the Mumford- 
Shah model seeks an optimal piecewise smooth function with smooth bound¬ 
aries to approximate the observed image. 

However, the Mumford-Shah model is hard to implement in practice be¬ 
cause the discretization of the unknown set of boundaries is very complex. 
Therefore, a parametric/explicit active contour method is used to represent the 
segmentation boundaries m • In addition, the special Mumford-Shah model 
with a piecewise constant approximation is studied by Chan and Vese m, and 
the level set method [371 is applied to solve this problem. Using an implicit 
representation of boundaries, the level set method has many advantages over 
methods using explicit representations of boundaries. For instance, the level 
set method handles the topological change of zero level set automatically [3] 
nog. Both the parametric/explicit active contour method and the level set 
method assume that each pixel belongs to a unique region. An alternative 
way to represent various regions is to use fuzzy membership functions [33 
EM5], which describe the levels of possible membership. Fuzzy membership 
functions assume that each pixel can be in several regions simultaneously with 
probability in [0,1]. These probabilities satisfy two constraints: (i) nonnegativ¬ 
ity constraint, i.e., the membership functions are nonnegative at all pixels; (ii) 
sum-to-one constraint, i.e., the sum of all the membership functions at each 
pixel equals one. Then the length of boundaries can be approximated by the 
Total Variation (TV) of fuzzy membership functions. The main advantages of 
using fuzzy membership functions over other methods include: i) the energy 
functional is convex with respect to fuzzy membership functions, guarantee¬ 
ing the convergence and stability of many numerical optimization methods, ii) 
fuzzy membership function has a larger feasible set, and it is possible to find 
better segmentation results. 

For two-phase segmentation, where there are only two regions, we only need 
one level set function or one fuzzy membership function. Multiphase segmenta¬ 
tion is more challenging than two-phase segmentation since more variables and 
constraints are involved in representing multiple regions and their boundaries 
effectively. The two-phase Chan-Vese model m has been generalized to multi¬ 
phase segmentation by using multiple level set functions to represent multiple 
regions [33]. Partitioning an image into N disjoint regions requires log 2 N level 
set functions. The advantage of using multiple level set functions is that it au¬ 
tomatically avoids the problems of vacuum and overlap of regions. However, 
the implementation is not easy, and special numerical schemes are needed to 
ensure stability [321l33U42U27j . For fuzzy membership functions, the sum-to-one 
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constraint is not satisfied automatically in multiphase segmentation. However, 
this constraint is easy to deal with in many cases, e.g., Fuzzy C-Mean (FCM) 
and its variants have closed-form solutions for the membership functions and 
are widely used in data mining and medical image segmentation [3 | [T | l38 p 6l 
129113 ll!30lf33[fT] . Variable splitting schemes are used in both [30] and [3T] to 
get efficient numerical algorithms. The Alternating Direction Method of Mul¬ 
tipliers (ADMM) method is used in [21] to derive the algorithm with two sets 
of extra variables. Primal-dual type algorithm is derived in m to solve the 
TV regularized FCM segmentation model. Both [21] and mu use projection 
to simple to handle the constraints of membership functions. Other segmen¬ 
tation approaches include a convex approach [9], two-stage methods mm, 
one single level set function approach m, et.al. 

Noise is unavoidable in images, and it is important to develop segmentation 
methods that work on noisy images. Among many types of noise, the Gaus¬ 
sian white noise is frequently assumed, and the L2-norm fidelity is adopted. 
However, when images are corrupted by non-Gaussian noise, in order to obtain 
a faithful segmentation, one has to model the noise according to its specific 
type |40U36U12l[23 yT 71H8lf7] . Particularly, the L 1-norm fidelity is used for both 
salt-and-pepper impulse noise and random-valued impulse noise in image pro¬ 
cessing mmm- In addition, it is robust to outliers and able to preserve 
contrast because the denoising process of L 1-norm models is determined by 
the geometry such as area and length rather than the contrast in the L2-norm 
case [T3] . 

Inspired by the fact that L 1-norm is more robust to impulse noise and 
outliers and can better preserve contrast, in this paper, we propose a vari¬ 
ational multiphase fuzzy segmentation model based on L 1-norm fidelity and 
fuzzy membership functions. This model can also deal with missing data in 
images. When there are missing pixels in an image, we randomly assign 0 or 
255 at these pixels by considering these pixels as corrupted by salt-and-pepper 
impulse noise. ADMM E2QIIX2I, which was rediscovered as split Bregman [21] 
and found to be very useful for LI and TV type optimization problems, is 
applied to solve this nonconvex optimization problem. By introducing two 
sets of auxiliary variables, we derive an efficient algorithm with all the sub¬ 
problems having closed-form solutions. In the theoretical aspect, we prove the 
existence of the minimizer and analyze the convergence of the algorithm. We 
note that the proposed method is closely related to the method in m since 
both methods use TV regularization and ADMM algorithm. The difference is 
that L 1-norm fidelity is considered in the proposed method, while L2-norm 
fidelity is used in [21]. 

The outline of this paper is as follows. In Section [2] we review some closely 
related existing works. In Section [31 the proposed model is described in detail, 
and the existence of the minimizer to the model is proved. In Section d] a 
numerical algorithm is derived, and its convergence analysis is presented. In 
Section experimental results and comparisons are presented to illustrate 
the effectiveness of the proposed method. Finally, the paper is concluded in 
Section EJ 
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2 Related works 


Let Q C M 2 be a bounded open subset with Lipschitz boundary, and let 
I : f2 be the given clean or noisy image. Let s = 1 for grayscale images 

and s = 3 for color images. Our goal is to find an 7V-phase “optimal” partition 
{^i}iLi such that Qi P| fij = 0 for all i ^ j and (J i=1 Q Define the set 
of TV-phase fuzzy membership functions as 


N 

A:= \ (ui, ..., un ) ui G BV(Q), Ui{x) > 0, ^^Ui(x) = 1, Vx G i? 

i=l 


where BV(fi) is the space of functions with bounded variation [2|. The closely 
related works are listed in the following and will be compared with our pro¬ 
posed method in Section [5J For the sake of simplicity, we use the notations 
U = (ui, • • • , un) and C = (ci, • • • , c/v) G M sAr , where q G M for grayscale 
images and q G M 3 for color images. 


— FCM [3]- Fuzzy c-means clustering method that solves 


min 

(' U,C)£AxR sN 



| I(x) — Ci \ 2 u?{x) dx 


using the alternating minimization algorithm. Though p can be any number 
no smaller than 1, it is commonly set to 2. 

— FCM_S2 [TB] - A variant of FCM that solves 


mm 

(U,C)G^xI 


N 

£ 

i =1 


Q 


| I{x) — Ci\ 2 u?(x) dx + a 


! I/(*) 

Js1 


°(x) dx 


where / is obtained by applying the median filter on I and a > 0 is a 
weight parameter. It can also be solved by the alternating minimization 
algorithm, and it is more robust to impulse noise than FCM. 

— FLICM |29j - Fuzzy Local Information C-Means clustering method that 
solves 


min 

(U,C )£AxR sN 



I I{x) — Ci | 2 v%(x) dx 



where A f(x) is a neighborhood of x. FLICM is more robust to both Gaus¬ 
sian noise and impulse noise than FCM. 

— L2FS [23j - L2-norm fidelity based Fuzzy Segmentation method that solves 


mm 

(' U,C)eAxR sN 


N 

£ 


:{/ \\Vui(x)\\dx + \ J \I(x) - Ci \ 2 Ui(x) dx^ 


( 1 ) 
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by ADMM. Here A > 0 is a parameter and f Q ||Viq(:r)|| dx denotes the TV 
of Ui with \\\7ui(x)\\ := y/(\7 Xl Ui(x)) 2 + ( \7 X2 Ui(x )) 2 for x = (x 1 ,x 2 ). For 
fixed C, [23] is related to the popular TV denoising method [39]. Note that 
the similar model is solved by other fast numerical methods in [3D]. 

— L1SS [26] - L 1-norm fidelity based Soft Segmentation method, in which 
log 2 N soft functions are introduced to represent N phases. Since the model 
for multiphase segmentation is complicated for more than four phases, we 
show the four-phase model as follows: 


mm 

ui,u 2 e[o,i],ceM siV 



HVuiC 


c)iicfe+Ay] f i 

i=l Jn 


| I{x) — Cj | Mj {x)dx 


where the membership functions Mj, j = 1, • • • ,4, are represented by soft 
functions ui(x),u 2 (x) G [0,1] in the following way: 


Mi(x) = ui(x)u 2 (x), M 2 {x) = ui(x)(l - u 2 (x)), 

M s (x) = (1 - u 1 (x))u 2 (x), M a (x) = (1 - ^i(x))(l - u 2 (x)). 


- L2L0 ^ - L2-norm fidelity and LO-norm regularization based image par¬ 
tition model: 

min ||Vm||o + X\\u - 1 1| 2 . 

u 

This model seeks a piecewise constant approximation of the original image 
/. Since this model can not specify the number of classes, a second step 
is applied to combine some classes if this model returns more classes than 
required. Here we apply FCM on the piecewise constant approximation u 
to obtain the final segmentation result. 


Remark: There are two advantages of our proposed method over L1SS. Firstly, 
we use fuzzy membership functions to represent regions, where N fuzzy mem¬ 
bership functions are required for an V-phase segmentation. Hence, the solu¬ 
tion space is much larger than L1SS, which ensures the higher possibility to 
obtain optimal segmentation. Secondly, the proposed method can take use of 
other commonly used segmentation methods such as FCM to gain good initial¬ 
ization of fuzzy membership functions. Multiphase segmentation is sensitive 
to initialization, and a good initialization is very important for a successful 
segmentation. However, it is hard to use the existing segmentation methods 
to get a good initialization for soft membership functions in L1SS. 


3 The proposed model 

In this paper, we assume that the given image can be approximated by a 
piecewise constant function, i.e., 


N 

7 ( x ) = +n(x). 

i=1 
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Here \ denotes the indicator function of region fii, i.e., 

( N / if xeOi\ 

Xttiyz) | q, otherwise, 

Ci is a constant that represents the given data in region and n can be 
outliers, impulse noise or other mixture types rather than Gaussian noise. 

Instead of using the L2-norm fidelity to measure the approximation error 
when the noise is the Gaussian white noise, we use the L 1-norm fidelity. Same 
as in the Mumford-Shah model, we require the segmentation boundaries to be 
smooth. Then we have the following model 


min 


E [ \\Vxsi i (x)\\dx + A f 
i—l J ^ J ^ 


N 




i= 1 


CiXni(x) 


dx , 


( 2 ) 


where A > 0 is a tuning parameter. Note that the TV of xot in the first term 
is equal to the length of boundary dQi. An equivalent form of model © is 


.mi 1 ! V( / W v Xni(x)\\dx + A f \I(x) - c*| xn^dxX . (3) 

{Hi}, c. =i \Jn Jq J 

Because \Q i can take values 0 and 1 only and {xf2i} is a partition, (xq 1 ,..., x^n) 
belongs to the set 

r N 

A := l (ui, ...,u N )\ui e BV{Q),Ui(x) e {0,1},^ Ui(x) = l,Vx E f2 

{ i=i 

At any point x E i?, there is only one function having value 1, and all the 
other functions have value 0. Thus set Aq is not continuous, which leads to 
difficulties and instabilities in numerical implementations. However, we can 
relax binary indicator functions {xnA to fuzzy membership functions {iq}, 
which satisfy the nonnegativity constraint and the sum-to-one constraint, i.e., 
(iq, ...,un) belongs to the set A defined in ([2|). It is obvious that Ui(x) E [0,1] 
and A is a simplex at any x E Q. So Ui(x) can be understood as the probability 
of pixel x to be in the ith class. Then we can rewrite our model © as 


min E(U,C)=H / II'Viz* (a) II dx + A f \I(x) - cAuAx) dxX . 

(u,c )eAxr* n Jn J 

(4) 

Note that model (pH) is convex with respect to U and C separately but not 
jointly. The difference between (pf]) and (pQ) is that the L2 fidelity term in (pQ) 
is replaced by the LI fidelity term. The existence of a minimizer for E(U, C) 
in (0) is stated in Theorem [lj 


Theorem 1 (Existence of a minimizer) Given an image I E there 

exists a minimizer ofE( U,C) in A x W N for any parameter A > 0. 
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Proof Since E(XJ, C) is positive and proper, the infhnum of E(XJ, C) must be 
finite. Let us pick a minimizing sequence (U n ,C n ) E A x M sAr that satisfies 
E(U n , C n ) infu,c E( U, C) as n oo. Then there exists a constant M > 0 
such that 


E(U n ,C n ) = Yj{J l|Vtt"(a:)|| dx + X J \I(x) - c?| <(®) dx J < M. 

Then each term in E( U n , C n ) is bounded, i.e., for each im 1, • • • , TV, 

f \\'Vu 7 l(x)\\dx < M, f 1 1 {pc) — Ci\u™{x)dx < M. (5) 

Jn JQ 


Since uf(x) E [0,1], we have f Q uf(x)dx < |i?|, where |i?| is the area of Q. 
Together with the first equality in (0, we have that uf is uniformly bounded 
in BV(fi) for all i — 1 ,..., N. By the compactness property of BV(f2) and the 
relative compactness of BV(f2) in L 1 (i?), up to a subsequence also denoted 
by {ui} n after relabeling, there exists a function u* E BV(Q) such that 

u™ — )> u* strongly in L 1 (i?), 
uf u* a.e. x E J2, 

Vu™ Vu* in the sense of measure. 


Then by the lower semicontinuity of the TV semi-norm, 


[ \\\7u*(x)\\dx < liminf f \\Vu^(x)\\dx. 
Jsi n ^°° Jn 


(6) 


Meanwhile since U n = (V™, ..., u 1 fo) E A, we have U* — (u*, ...,7qjy) E A. 

It is easy to derive the first order optimality condition about c™, which is 


0E / d\I(x) — c™\u™(x)dx, 

Jn 


where d\ • | is the sub differential of | • |. Since uf(x) > 0 and f Q uf(x)dx > 0, 
the above equation implies that the constant Pf satisfies 


\c?\ < 


— Ir lloo- 


By the boundedness of sequence {c- 1 }, we can extract a subsequence also de¬ 
noted by {Pf} such that for some constant c* 


c, , as n oo. 


Finally, since uf(x) —> u*(x), a.e. x E Q and Pf —>• c* as n oo, Fatou’s 
lemma yields 

f | I(x) — c*\u*(x)dx < liminf [ \I(x) — Pf\u 7 f{x)dx. (7) 

Jii n ^°° Jii 
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Combining inequalities © and ([7j) for all i, on a suitable subsequence, we 
established 


£(U*,C*) < liminf E(TJ n , C n ) = inf E( U,C), 

n—>■ oo U.C 

and hence (U*, C*) must be a minimizer of the energy E. This completes the 
proof. □ 

The minimizer of £'(U, C) is not unique due to the following hidden sym¬ 
metry property. Denote Sn as the permutation group of {1, ...,7V}, i.e., each 
permutation a G Sn is defined as a one-to-one map a : {1,..., TV} — >> {1,..., TV} 
such that {(j(l),..., <t(TV)} is a rearrangement of {1,...,TV}. Denote U a = 
(^o-(i), • • • ? ^cr(N )) 1 = (c a ( 1 ),..., c a (jv))- It is straightforward to show the fol¬ 
lowing theorem. 

Theorem 2 (Symmetry of minimizer) For any (U, C) G A x R N and 

any a G Sn, 

E(\J a , C a ) = E(XJ, C). 

In particular, suppose that (U*,C*) is a minimizer of E( U, C), i.e., 

(U*,C*)= argmin £(U,C). 

(U,C)eAxR N 

Then, for any a G Sn, we have 

(U;,C;)= argmin E( U,C). 

(u.cje^xK^ 


4 The numerical algorithm and its convergence analysis 

In this section, we provide an efficient algorithm based on ADMM and discuss 
its convergence. 


4.1 The algorithm 


ADMM is applied, in this subsection, to solve the proposed fuzzy segmentation 
model We introduce two sets of auxiliary variables D = (di,..., d/v), W = 
(wi, ...,wn) such that Vrq = d^,iq = Wi. Then the model ((TJ) is equivalent to 
the following minimization problem with equality constraints: 

D min u ^ \\di(x)\\dx + X J^\I(x) - c*| Wi(x) dx j + ^(W), ^ 

subject to Viq = di,Ui = Wi,\/i = 1,..., TV, 


where 5 a is the indicator function of the set A, i.e., 


(W) = | 


0, if W € A, 

+ 00 , otherwise. 
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The augmented Lagrangian function for problem ((HJ) is: 


if(D,W,C,U;A D ,A w ) 




JS 

= £< 
i =1 

yj \\di(x)\\ dx + X j 

f IFd- c il 

n 

Mx) il . 

|+^(W) 

N 1 

{Xdi 5 V 

-*)+i 

f |Vw,(x) 

JQ 

-di(x)f 

dx | 

N 1 

f 


r 

'I 


1=1 

(A wi •> ^i 

W ’ ) + r 2j 

/ 1 Ui(x) - u 
'n 

h(x )| dx 

> 5 


where tId = (A<q,A^), Hw = (A Wl , X WN ) are the Lagrangian multipliers 
and r is a positive constant. Here (A Vrq— = f Q A ^.(x)(Vui(x)—di(x))dx, 
and (X Wi , Ui - Wi) = f Q \ w . (x)(m(x) -Wi(x))dx. 

The ADMM solves the primal variables one by one in the Gauss-Seidel 
style and then updates the dual variables (Lagrangian multipliers). It can be 
summarized as follows: 

D /c+1 = argmin Jgf (D, W*, C k , U k ; A^), 

D 

W fe+1 = argminJ^(D fc+ \ W, C fc , 

w 

C fc+1 = argmin^(D fc+1 ,W /c+1 ,C,U /c ;yl^,yl^), 

c 

U^ 1 = arg min ££ (D /c+1 , W /c+1 , C /c+1 , U; A^), 

In the following, we show how to solve each subproblem and then give the 
algorithm. 


D-subproblem 

The subproblem for D is equivalent to 

D fc+1 = argminW < f \\di(x)\\dx + ^ [ 
D i=1 Ji? z Jn 


di{x) - - 


x dM) 


dx 


This is separable and the optimal solution of d k+1 can be explicitly expressed 
using shrinkage operators. We simply compute 


d*+\x)=S V«?(*) + 
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where S is the shrinkage operator defined as 

S(v,t) := ijL- *max(|M| -r,0). 
For the sake of simplicity, we denote this step as 

r r J 


D fc+1 = S ^VU fe 


W-subproblem 

The subproblem for W is equivalent to 

N 


mm 

w 


WL 


, . k, s X w( x ) X \I{x)-C k i\ , . - 

Wi{x) - u k (x) - ^ dx } + S A ( W). 


Since A is a convex simplex at any x G i?, the solution is given by 

N 


w fc+1 = n A 


x t, A|/-cf| 


i=lj 


where 11a is the projection onto the simplex A , for which more details can be 
found in m- We denote the step as 

w fe+1 = n A fu fe + 

r r 


C-subproblem 

The subproblem for C is 

TV 

C /c+1 = argminV^ / \I(x) — C{\w^ 1 (x) dx. 
c ^Jn 

It is separable, and c\ +1 can be solved independently. The optimality condition 
for each c^ +1 is 

0 G - [ d\I{x) - <$ +1 \wt +1 (x)dx. (9) 

JQ 

Because the right-hand side of (0 is maximal monotone [3], the bisection 
method and ADMM are applied to solve it EM- The next lemma provides 
a new way to find an optimal solution for discrete cases. 
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Lemma 1 Given a finite non-decreasing sequence i.e., 

I[ 1] < ^[2] < ••• < ^[n]j 

and a non-negative sequence {w^}f =1 with A = w [j] > 0; the optimal 

solution set for 

n 

™ n SlL] ~ c \ w iib ( 10 ) 

3= 1 

is [l\j*],I\j*+i\]\, where j* and j* satisfy 

r r -1 

A - 2j2 w [j] < 0 < A - 2 w [j]> 

3 =1 3=1 

3* + 1 3* 

A-2j2 w [j]<0< A ~ 2 J2 w bV 

3=1 3=1 

The fuzzy median of {Ij}f =1 with respect to the weight {w^p}f =1 J28 1/ . which 
is defined as (I[j*] + %*+i]])/2, is an optimal solution. If, in addition, there 
exists j * such that 


3 * 3 *~ 1 

A - 2 M, b1 <0<A - 2 H w [i]> 

J = 1 J = 1 

then j* = j* — 1, and the optimal solution is unique. 

Proof The optimality condition of m is 

n 

0 € h(c) := J2 d \ I W ~ c \ w bV 

3 = 1 

We can see that h(c) is non-increasing and it can be expressed as 


(A, 

[A - 2w m ,A\, 
A - 2w[q, 


c < -f[i] ? 


c = I\ 


[i]» 


c e (/[i],/[ 2 ]), 


h(c) = < 


[A-2 i wy], A - 2 ELi w [i]] > c = J [s] > 

A ~ 2 Ylj=i w ij}, (k e (i[ 8 ],i[ s+ 1 ]), 


-A 


c > L n ]. 


The range of h is [—A, A], and h(c) can take multiple values when c = Jyj for 
any j with wy] ^ 0. Therefore, we can always find j* such that 

3* 3*~ 1 

A — 2 W[j] <0<A-2^ W[i] . 

J=1 j=l 
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Thus 0 G /l(/[j*]) = 


[- 4-2 Y$=i w m > A ~ 2 Ejh 1 w\ji 

mal solution. In addition, we have that h(c) >0 when c < I[j 
can find j* such that 


and /[j*] is an opti- 
Similarly, we 


3* 

A ~ 2 Y1 w \j\ < 0^-4 -2 J2 w li}’ 

3=1 3=1 

and I[j*+ 1 ] is an optimal solution. In addition, we have that h(c) < 0 when 
c > /[j # -|4]. Then h(c) being non-increasing gives us that the set of optimal 
solutions for (hud is 7[j # _|_ij]. When j* = j* + 1, the optimal solution is 

unique. □ 

Remark: When there are missing pixels in images, we can put a mask on 
the data fidelity term as in image inpainting problems [31] or assign a value 
to each missing pixel. In [23], the missing pixels are assigned with zero, and C 
changes based on the percentage of missing pixels. While, this lemma tells us 
that assigning 0 or 255 (for grayscale images) randomly to these missing pixels 
will not change the value of c with a high probability because q is a weighted 
median. Also, by assigning 0 or 255 randomly, this algorithm is able to deal 
with cases where more than half of the pixels are missing. See the numerical 
experiments in Section 0 

Assume that I and w^ 1 are vectors in M n , where n is the total number 
of pixels. We can reorder the components of / and w^ 1 such that the mono¬ 
tonicity of / in Lemma |T] is satisfied. If there are multiple optimal solutions 
for q, we choose the smallest one. We denote this step as 

C^ 1 = ^(W fe+1 ). 


U-subproblem 

The subproblem for U is equivalent to 

N 


U fe+1 = argmin^^ [ 
u “f 


i= 1 


Q 


Vui(x) - dj +1 (x ) + 


X i( x ) 


dx 


s 

J Q 


Ui{x) — u^ +1 (x) + 




dx. 


It is separable for r^ +1 , and, from the following first order optimality condition 
for each u^ 1 : 


V T V^ + 1 (x)-^ + 1 (x) + 


AM) + („*+!(*) - W *+\ X ) + MM 


= o, 
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we can derive the closed-form solution of u ^ +1 from the equation: 

(V T V + I) u k+1 (x) = V T d k+1 (x) + w k+1 (x) - - - AiM. 

A diagonalization trick can be applied to find u ^ +1 efficiently [35] , We denote 
this step as 

u fe+1 = (V T V + I) -1 fv r D fc+1 + w fe+1 - _ dw 

v / y r r 

Updating dual variables. We denote the steps as 

A'g 1 =A k 10> + r (VU fe+1 - D fe+1 ), 

A** 1 =A k w + r (U fc+1 - W fe+1 ). 

Combining all these steps together, we obtain the proposed LI Fuzzy Seg¬ 
mentation algorithm (L1FS) in Algorithm 1 for solving (|2j). 


Algorithm 1 The proposed L1FS algorithm 

— Initialization: U° and C° are specified, = 0 , = 0 . 

— For fc = 0,1, 2, • • •, repeat until the stopping criterion is reached. 


D /c+1 

W ^ 1 

C k+i 

U fc+1 

A k n +1 


s(v\J k + A±, - 

\ r r 


n A [u k + 


A ^ A|/-C fc | 


ip(W k+1 ), 


(V T V + I) 1 ( V r D 


vT-n/c+i 


W 


■fc+i 


A^ + r(VU k+1 -D k+1 ), 

A^ + r(U fe+1 - W fc+1 ). 





- Output: C fc+1 ,U fc+1 . 


Remark: Though there are four variables D, W, C and U, they can be 
grouped into two variables (D, W) and (C,U) and the subproblems for these 
two variables can be decoupled into the four subproblems. So it is essentially 
a two block ADMM applied on a nonconvex optimization problem. 

Because problem is nonconvex, a good initial guess of U° and C°, 
which can be obtained from other methods using fuzzy membership functions, 
is helpful for the stability of L1FS. Thus we update D 1 and W 1 first because 
both of them can use the initial guess of U°. 
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4.2 Convergence analysis 

If C is given, problem (|2|) is convex. Then, the algorithm is a standard ADMM 
by considering (D, W) together, and its convergence is guaranteed [5]. In this 
section, we give a convergence result of Algorithm 1 for unknown C. To sim¬ 
plify notations, let us define the sextuples: 

Z := (D,W,C,U,A d ,A w ). 

A point Z is a KKT point of ([H]) if it satisfies the following KKT conditions: 


d\\d*(x)\\-\* di (x)3 0, 

(Ha) 

X\I-C\ — A w + d5 A (W) 3 0 , 

(lib) 

f d\I{x) — c*\w*(x)dx 3 0, 

J Q 

(11c) 

K T X di (x ) + X Wi (x) = 0, 

(lid) 

\7ui(x) - di(x) = 0, 

(lie) 

Ui(x) — Wi(x) - 0. 

(Ilf) 


where <9|| • || and are sub differentials of || • || and 5^(-), respectively. 

Theorem 3 (Convergence to stationary points of L1FS) Let {Z /c }^ 1 
be a sequence generated by Algorithm 1 that satisfies the condition 

lim (Z k+1 - Z k ) = 0. 

k-^oo V 7 


Then any accumulation point of {Z /c }^ 1 is a KKT point of problem J3p. 
Consequently, any accumulation point of (G k ,XJ k ) is a stationary point of 
problem W- 


Proof From the updating formulas in Algorithm 1, we obtain the following 
inequalities related to the subsequent iteration: 


w k+1 -w k = n A (u k 


V T A$- 


+ 1e, 

r r J 

(12a) 

+ A / — C fc |\ w „ 

r r ) 

(12b) 

)~c k , 

(12c) 

1 - VU fc ) + (W fe+1 - U fe ) 

(12d) 

) 

r ’ 

(12e) 

1 _ D fc+1 ) , 

(12f) 

- W fc+1 ). 

(12g) 
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By the assumption lim (Z /c+1 — Z k ) = 0, the left-hand side and right- 

k—too 

hand side of each equality in m go to zero as k oo. Then we have 

5 (vv k + ^, ^ - D fe -> 0, (13a) 

n A ^U fe + ^ - - w k -> 0 , ( 13 b) 

V>(w fc+1 ) - C k -> 0, (13c) 

V r (D fc+1 - VU fc ) + (W fc+1 - U fc ) - ZMd _ 3* 0 , (i 3d ) 

r r 

VU fe+1 - D fe+1 -> 0, (13e) 

U fe+1 - W fe+1 -> 0. (13f) 


Assume Z* = ( D *, W *, C *, U *, A ^. A ) is an accumulation point of Z fc . 
m gives us 


n A u 


5 (vU* + - 

\ r r 

! AjJ - C 

r r 


-D* = 0, 
- W* = 0, 


^(w*) - c* = o, 

V T ^ + A^f = 0, 

VU* - D* = 0 5 

u* - W* = 0 . 

By (114a.ll . it follows that D* is a solution of the minimization problem 


(14a) 

(14b) 

(14c) 

(14d) 

(14e) 

(14f) 


r^t{fj d , M ld x+ r -f a 


di{x) - \7u*(x) 


KM) 


dx 


Thus D* satisfies the first order optimality condition 
0 E 3||d*(x)|| + r (d*(x) — \7u*(x) — 


KM) 


Using (lUel) . we simplify the above condition as 

o g < 9 |K(x)|| - KM)- 

By G20, W* is a solution of the following minimization problem: 


(15) 


min 

w 



Wi(x) - u*{x) - 


KM) , A|j(x)-c?| 



+ 5 a (W). 
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Hence W* satisfies the optimality condition 


0 G dS A ( W*) + r(W* - U*) - Aw + X\I - C*|. 

Together with the equality in (114fjh the above equation can be simplified as 


The equation 


0 G dS A { W*) - + X\I - C*|. (16) 

cj) implies that C* is a solution of the minimization problem 

N 


i n T / \i( x )- c i\ w *( x ) dx, 


mm 

c 


and the optimal condition is 

OG / d\I(x) — c*\w*(x)dx. (17) 

Jn 

The equivalence of equations (I14al) - (ll4cl) with (fTSl)-(H7l). together with 
equations (|14d|) - (ll4f|) shows that the accumulation point Z* satisfies the KKT 
condition in equations (lllal) - (lllflh thus Z* is a KKT point of problem 

Since problem © and problem © are equivalent, the convergence of 
(C^TJk) in the statement follows directly. □ 

The convergence analysis is motivated by uni- This theorem implies that 
whenever {Z /c }^ 1 converges, it converges to a KKT point of problem ©• 
However, since the minimization problem is nonconvex, the KKT point is 
not necessary to be an optimal solution. 


5 Experimental results 

In order to demonstrate the effectiveness of the proposed method, we compare 
our method with some existing methods on both synthetic and real images. 
These methods (FCM, FCM_S2, FLICM, L2FS, L2L0, and L1SS) are discussed 
in Section [21 Since the segmentation models in these methods are not jointly 
convex, and they are easy to get stuck in local minima, “good” initialization 
is crucial for all algorithms, especially when the given image is noisy. For 
FCM, FCM_S2, and FLICM, the initial U is uniformly distributed in [0,1] and 
normalized to satisfy the sum-to-one constraint. While for TV based methods 
L2FS and L1FS, one can also use the results of FCM and FCM_S2 as the initial 
U and C. Here we consider three ways for choosing the initial U and C: the 
result of FCM, the result of FCM_S2, and U as functions uniformly distributed 
in [0,1] and C as the result of FCM. In all the experiments, we choose the one 
with the highest performance among all the three initializations. For L1SS, we 
use the initialization method as described in the original paper. The stopping 
criterion, which is the same for all the methods except L1SS, is defined as 

||U fe+1 — U fc || 2 ^ 

Ilu*ll2 <e 
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where e is a very small number. For L1SS, this stopping criterion does not work 
since L1SS leads to contrast loss due to the error in calculating class centers 
{ci}f =1 at early iterations. However, the contrast of L1SS will be enhanced 
gradually if the number of iterations increases. To gain satisfactory results, we 
choose to stop L1SS by setting the maximum number of iterations to be 1000. 

The parameters of the methods being compared are set as follows. In FCM, 
we set p = 2. In FCM_S2, we set p = 2, a = 5, and the window size for the 
median filter as 5 x 5. In FLICM, we set p = 2 and the window size of the 
neighborhoods as 3 x 3 or 5 x 5, which depends on the noise level. However, the 
weight parameter A for L2FS, L1SS, and L1FS are tuned for each experiment 
to achieve optimal results. In all experiments, the range of A for L2FS is 
[0.00005,0.0005], for L1SS is [0.03,1], and for L1FS is [0.001,0.05]. For all 
methods, e = 10 -6 for the two-phase segmentation and e = 10 -4 for the 
multiphase segmentation. We use the default parameters in L2L0 except that 
we set A = 10. 

The clustering results of all methods are obtained by checking the max¬ 
imum value of their membership functions. Then we display the recovered 
piecewise constant image as the final result. To compare segmentation results 
quantitatively, we consider the Segmentation Accuracy (SA) defined by 

SA # correctl y classified pixels 
#all pixels 



(a) (b) (c) 


Fig. 1 Synthetic piecewise constant test images, (a) A two-phase grayscale image with 
intensities 20 and 128, size 128 X 128; (b) A five-phase grayscale image with intensities 0, 
63, 127, 192, and 255, size 235 x 237; (c) A six-phase color image with color vectors (12 11 
242), (242 12 11), (242 241 242), (243 241 12), (12 12 12), and (12 242 12), size 100 X 100. 

The synthetic piecewise constant test images are displayed in Fig. HJ The 
noisy images are contaminated by three types of noise: Gaussian Noise (GN) 
with zero mean and standard deviation varying from 10 to 80, Salt-and-Pepper 
Impulse Noise (SPIN) with 10% to 60% pixels contaminated, and Random- 
Valued Impulse Noise (RVIN) with 10% and 60% pixels contaminated. 


5.1 Test on Fig. flal 

Tab. [U provides the SA comparison of these six algorithms, and Figs. [20 show 
some of the corresponding results. 
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Table 1 SA performance of different methods applied on Fig. Ilal contaminated by dif¬ 
ferent levels of GN, SPIN, and RVIN. 


GN (cr) 

10 

20 

30 

40 

50 

60 

70 

80 

FCM 

1 

0.9970 

0.9642 

0.9146 

0.8627 

0.8162 

0.7811 

0.7524 

FLICM 

1 

0.9999 

0.9996 

0.9990 

0.9975 

0.9968 

0.9960 

0.9954 

L2FS 

1 

1 

1 

1 

1 

1 

0.9999 

0.9998 

L1SS 

1 

1 

1 

0.9998 

0.9996 

0.9984 

0.9975 

0.9965 

L1FS 

1 

1 

1 

1 

1 

1 

1 

0.9998 

SPIN (%) 

10% 

20% 

30% 

40% 

50% 

60% 



FCM 

0.9480 

0.8983 

0.8478 

0.7979 

0.7486 

0.6980 



FLICM 

0.9984 

0.9921 

0.9738 

0.8002 

0.7313 

0.6554 



L2FS 

0.9999 

0.9998 

0.9982 

0.9983 

- 

- 



L1SS 

0.9998 

0.9990 

0.9977 

0.9967 

0.9956 

0.9953 



L1FS 

1 

1 

1 

1 

1 

0.9995 



RVIN (%) 

10% 

20% 

30% 

40% 

50% 

60% 



FCM.S2 

0.9985 

0.9972 

0.9945 

0.9862 

0.9630 

0.9042 



FLICM 

0.9987 

0.9985 

0.9970 

0.9958 

0.9948 

0.9919 



L2FS 

1 

1 

1 

0.9998 

0.9995 

0.9974 



L1SS 

1 

0.9995 

0.9985 

0.9979 

0.9966 

0.9891 



L1FS 

1 

1 

1 

1 

1 

0.9976 





(a) GN a=30 (b) 0.9642 (c) 0.9996 (d) 1 (e) 1 (f) 1 



Fig. 2 Two-phase segmentation on the synthetic image Fig. Ilal with Gaussian noise. First 
column: images contaminated by Gaussian noise with standard deviations 30 and 60, re¬ 
spectively; Second column to last column: results of FCM, FLICM, L2FS, L1SS, and L1FS, 
respectively. The SA values are reported below each segmentation result. 


For GN, all methods being tested give correct segmentation results for 
small standard deviations (e.g., a = 10). As the standard deviation increases, 
the SA value of FCM decays very fast. All the other algorithms have very 
large SA values even when the standard deviation is 80. L1FS has the best 
performance for all cases, and it is able to give correct segmentation results 
even when a < 70. L2FS is the second best algorithm, and it is able to give 
correct segmentation results when a < 60. 

The results of all methods when a = 30 and 60 are displayed in Fig. [21 The 
results of FCM are relatively “noisy” (the second column). For FLICM (the 
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(a) SPIN 20% (b) 0.8983 (c) 0.9921 (d) 0.9998 (e) 0.9990 (f) 1 



(g) SPIN 40% (h) 0.7979 (i) 0.8002 (j) 0.9983 (k) 0.9967 (1) 1 


Fig. 3 Two-phase segmentation on the synthetic image Fig. Ilal with SPIN. First column: 
images contaminated by 20% and 40% SPIN, respectively; Second column to last column: 
results of FCM, FLICM, L2FS, L1SS, and L1FS, respectively. The SA values are reported 
below each segmentation result. 




Fig. 4 Two-phase segmentation on the synthetic image Fig. Ilal with RVIN. First column: 
images contaminated by 20% and 40% RVIN, respectively; Second column to last column: 
results of FCM_S2, FLICM, L2FS, L1SS, and L1FS, respectively. The SA values are reported 
below each segmentation result. 


third column) and L1SS (the fifth column), the segmentation error occurs on 
the middle edge. 

For SPIN, Tab. [1] shows that FLICM performs much better than FCM for 
noise levels 10%-30%. However, if the noise level is higher than 30%, both 
FCM and FLICM have very poor performance. L1SS achieves much better 
performance than FLICM, even when the noise level is higher than 30%. L2FS 
is slightly better than L1SS. Meanwhile, if the noise level is higher than 50%, 
L1SS fails to give a satisfactory result. L1FS has the highest SA among all 
methods. It gives completely correct segmentation results for noise levels 10%- 
50%. Fig. 0 shows the results of all methods for noise levels 20% and 40%, 
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respectively. For L2FS and L1SS, the segmentation errors occur along the 
middle edge. We also observe that for high noise levels such as 40%, both L2FS 
and L1SS suffer from slight contrast loss, e.g., Fig. |3j] and Fig. [30 However, 
L1FS is still able to preserve contrast, e.g., Fig. [33 and Fig. [511 

For RVIN, FCM_S2 is used to initialize U and C for TV based methods 
L2FS and L1FS. Tab.[T]shows that FCM_S2 has the worst performance among 
all. L1SS is much better than FCM_S2 especially for high noise levels. FLICM 
performs slightly better than FCM_S2. L1FS achieves the best performance 
which is slightly better than L2FS. L2FS gives correct segmentation results at 
noise levels 10%-30%, while our method L1FS can give the correct segmenta¬ 
tion results at noise levels 10%-50%. Fig. 0] shows the results for noise levels 
20% and 40%. Again, we find that most of the errors occur along the middle 
edge for FLICM and L1SS. Moreover, the results of L2FS in Fig. [TJ] L1SS in 
Fig. 0e] and Fig. 0k] lose some contrast. 

Next we analyze the contrast problem for TV based methods. For L2FS, 
the estimated intensity Ci in each segmented region roughly equals the mean 
value of the intensities in that region. In the Gaussian noise case, the noise 
has zero mean and therefore it has almost no impact on the estimation of q. 
However, for both impulse noise cases, the noise has a significant influence 
on the estimation of q by taking the average. More specifically, assuming 
that the impulse noise follows the uniform distribution, its impact on the 
estimation of c* is like this. Given an image with intensity range [0,255], for 
the region Qi with true intensity /3 < 128, if there are more noisy pixels with 
intensity greater than /3 than those with intensity less than /?, then q > /? 
after taking the average and vice versa. Hence, the range of the image will be 
shrunk by applying L2FS even when the segmentation is correct, and thereby 
the recovered image will suffer from contrast loss. Note that the contrast loss 
problem has also been reported for the TVL2 restoration model in impulse 
noise removal [f2j. For L1SS algorithm, one step of ADMM is used to solve 
the C-subproblem approximately. Thus q is not accurate for the first few 
iterations. However, since L1SS uses the L 1-norm fidelity, the loss of contrast 
becomes more and more subtle as the number of iterations increases. In L1FS, 
we solve the C-subproblem exactly. Thus, L1FS can preserve contrast well in 
the segmentation process. 


5.2 Test on Fig. [lb] 

The performance comparison for the multiphase synthetic piecewise constant 
gray image Fig. [lb] is shown in Tab. [2] and Fig. 01 

As shown in Tab. 01 FCM performs poorly for GN, while L2FS and L1FS 
perform relatively well with similar SA. For the noise level a = 10, both L2FS 
and L1FS give a correct segmentation result. For SPIN, FCM also gives the 
worst performance in terms of SA, while L1FS achieves the best performance. 
L2FS fails to yield a correct segmentation result when noise levels a > 50. For 
RVIN, FCM_S2 achieves high SA values since it can smooth out some noise 
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Table 2 SA comparison of different methods applied on Fig. rm with different levels of GN, 
SPIN, and RVIN. 


GN (cr) 

10 

20 

30 

40 

50 

60 

70 

80 

FCM 

0.9987 

0.8191 

0.6634 

0.5849 

0.5233 

0.4718 

0.4319 

0.4017 

L2FS 

1 

0.9999 

0.9994 

0.9978 

0.9959 

0.9950 

0.9931 

0.9918 

L1FS 

1 

0.9999 

0.9993 

0.9980 

0.9964 

0.9950 

0.9931 

0.9905 

SPIN (%) 

10% 

20% 

30% 

40% 

50% 

60% 

- 

- 

FCM 

0.9202 

0.8431 

0.7638 

0.6847 

0.6096 

0.5296 

_ 

_ 

L2FS 

0.9926 

0.9877 

0.9713 

0.9673 

- 

- 

- 

- 

L1FS 

0.9977 

0.9948 

0.9923 

0.9894 

0.9848 

0.9782 

- 

- 

RVIN (%) 

10% 

20% 

30% 

40% 

- 

- 

- 

- 

FCM 

0.9922 

0.9809 

0.96672 

0.9248 

_ 

_ 

_ 

_ 

L2FS 

0.9949 

0.9923 

0.9880 

0.9731 

- 

- 

- 

- 

L1FS 

0.9976 

0.9957 

0.9922 

0.9868- 

- 

- 

- 



Table 3 SA comparison of different methods applied on Fig, flcl with different levels of GN, 
SPIN, and RVIN. 


GN (cr) 

10 

20 

30 

40 

50 

60 

70 

80 

FCM 

1 

1 

0.9998 

0.9958 

0.7927 

0.7772 

0.7488 

0.7205 

L2FS 

1 

1 

1 

1 

0.9996 

0.9992 

0.9989 

0.9978 

L2L0 

1 

1 

1 

0.9999 

0.9996 

0.9991 

0.9983 

0.9967 

L1FS 

1 

1 

1 

1 

0.9998 

0.9994 

0.9985 

0.9973 

SPIN (%) 

10% 

20% 

30% 

40% 

50% 

60% 

- 

- 

FCM 

0.8498 

0.7294 

0.6128 

0.5092 

0.4248 

0.3501 

_ 

_ 

L2FS 

0.9960 

0.9925 

0.9883 

0.9822 

0.9772 

- 

- 

- 

L2L0 

0.9951 

0.9880 

0.9819 

0.9740 

0.9401 

0.8752 

- 

- 

L1FS 

0.9973 

0.9937 

0.9897 

0.9854 

0.9810 

0.9732 

- 

- 

RVIN (%) 

10% 

20% 

30% 

40% 

50% 

60% 

- 

- 

FCM 

0.8971 

0.7992 

0.6967 

0.6085 

0.5196 

0.4294 

_ 

_ 

L2FS 

0.9974 

0.9957 

0.9899 

0.9853 

0.9841 

0.8988 

- 

- 

L2L0 

0.9971 

0.9955 

0.9906 

0.9856 

0.9727 

0.9488 

- 

- 

L1FS 

0.9988 

0.9963 

0.9939 

0.9906 

0.9881 

0.9777 

- 

- 


in the segmentation process. L2FS performs much better than FCM_S2, and 
L1FS performs the best among all methods. 

Fig.[5]shows some results corresponding to Tab|2l The first row is the noisy 
images being tested. The second row shows the results of FCM (Fig. |5gpj| ) 
and FCM_S2 (Fig. I5kll5l)h Most of them looks very “noisy” except Fig. (5kj 
L2FS and L1FS give very clean results in the third row and last row, respec¬ 
tively. Comparing the results by L2FS and L1FS for Gaussian noise, both 
of them have high visual qualities. However, for SPIN and RVIN, it is obvi¬ 
ous that L1FS preserves the contrast much better than L2FS and has better 
segmentation results. 
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(g) 0.6634 (h) 0.4718 (i) 0.8431 (j) 0.6847 (k) 0.9809 (1) 0.9248 




Fig. 5 Multiphase segmentation on the synthetic image Fig. Ilbl with different levels of GN, 
SPIN, and RVIN. First row: images contaminated different types of noise with different 
levels. Second row to last row: results of FCM, L2FS and L1FS, respectively. The SA values 
are reported below each segmentation result. 


5.3 Test on Fig. [Tel 

In Tab. [3] and Fig. El we test the multiphase synthetic piecewise constant color 
image Fig. [T^ with various levels of GN, SPIN, and RVIN. 

From Tab. El in the GN case, when the standard deviation of noise a < 20, 
all the four methods, including FCM, L2FS, L2L0 and L1FS, give correct seg¬ 
mentation results. Moreover, both L2FS and L1FS yield correct segmentation 
results when a < 40. When a > 50, the performance of FCM decreases rapidly, 
while L2FS, L2L0, and L1FS still achieve very large SA values. Note that we 
initialize U randomly for GN in this test. For SPIN and RVIN, FCM has the 
worst performance which is far lower than that of the other three methods. It 
is also obvious that L1FS outperforms L2L0 and L2FS. 

Fig. El shows some results corresponding to Tab. El The first row shows 
the tested noisy images. The results of FCM in the second row seems to be 
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IBV 


GN <r=30 

GN (7=60 SPIN 20% 

SPIN 40% RVIN 20% 

RVIN 40% 


f . 

wskti /.vi 
jj j 

A 


0.9998 0.7772 0.7294 0.5092 0.7992 0.6085 



1 0.9992 0.9925 0.9822 0.9957 0.9853 



1 0.9991 0.9880 0.9740 0.9955 0.9856 



1 0.9994 0.9937 0.9854 0.9963 0.9906 


Fig. 6 Multiphase segmentation on the synthetic color image Fig. llcl with different levels of 
GN, SPIN, and RVIN. First row: images contaminated different types of noise with different 
levels. Second row to last row: results of FCM, L2FS, L2L0, and L1FS, respectively. The SA 
values are reported below each segmentation result. 


“noisy” in most cases. The results of L2FS (the third row), L2L0 (the fourth 
row), and L1FS (the last row) are very clean. However, in terms of contrast, it 
is obvious that L1FS outperforms L2L0 and L2FS particularly for SPIN and 
RVIN. 
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5.4 Test on real images 

We test some real images including two medical images and six natural images 
in Fig. 0 without artificial noise. However, the images can be regarded as con¬ 
taining some types of noise due to the acquisition and transmission processes. 

The results of FCM and L1FS are displayed for comparison. One can see 
that FCM tends to produce some tiny components and irregular segmenta¬ 
tion boundaries. By contrast, L1FS tends to smooth out tiny components to 
generate large ones and smooth boundaries between regions, which is more 
natural for human vision system and good for postprocessing. This smoothing 
effect is mainly achieved by total variation regularization in the L1FS model. 
Moreover, L1FS preserves slightly better contrast in the piecewise constant ap¬ 
proximation than FCM, which is mainly achieved by the use of the L 1-norm 
fidelity. 

In Figs. EP FCM and L1FS give quite different segmentation results. Ob¬ 
viously, FCM fails to segment the blue color part of the clothes in the original 
image, while the proposed L1FS works well. To illustrate the difference of these 
two methods, we display, in Fig. [HI the corresponding six segmented regions 
of FCM and L1FS for the women image Fig. [7]] respectively. We find that the 
segmented regions of FCM in Fig. lHallHel are somehow “noisy”. In particular, the 
background lattice pattern is partitioned into five regions as shown in Fig. [Hal 
IHel Compared with FCM, the proposed L1FS gives quite clean segmented 
regions in the second row in Fig.[8] Especially, the background lattice pattern 
are classified into only two classes as shown in Fig. [Hg[l8hl We further compare 
the blue parts of the clothes corresponding to Fig. [8c] and Fig. [HU In Fig. [8cl 
some background pattern heavily affects the estimation of C, and therefore 
the color is not blue any more. However, in Fig. [HU since the background is 
clean, the correct color can be obtained. To sum up, Fig. [H] demonstrates that 
the proposed L1FS gives smoother segmentations than FCM. 


6 Conclusion 

This paper presents a novel piecewise constant image segmentation model 
based on fuzzy membership functions and L 1-norm fidelity. ADMM is applied 
to derive an efficient numerical algorithm, in which each subproblem has a 
closed-form solution. In particular, an efficient algorithm is proposed to find 
the fuzzy median. The numerical results on both synthetic and real piecewise 
constant images demonstrate that the proposed method is superior to some 
existing state-of-the-art methods since it is more robust to impulse noise and 
can preserve better contrast. Even in the case of Gaussian noise, the proposed 
method can achieve similar results as its L2-fidelity counterpart. In this work, 
we assume that the image to be dealt with is piecewise constant, which works 
well on images with homogeneous image features. The future work is to extend 
this framework to piecewise smooth image segmentation. 
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Fig. 7 Segmentation on real images. First column and fourth column: real color images; 
Second column and fifth column: results of FCM; Third column and last column: results of 
L1FS. 



Fig. 8 Detailed comparison of FCM and L1FS on the woman image. First row: six seg¬ 
mented regions of FCM. Second row: six segmented regions of L1FS. 
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